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Abstract 

We describe a method to measure the concentration of 224 Ra and 226 Ra in the heavy water target used to detect solar 
neutrinos at the Sudbury Neutrino Observatory (SNO) and in the surrounding light water shielding. A water volume of (50- 
400) m 3 from the detector is passed through columns which contain beads coated with a compound of manganese oxide 
onto which the Ra dissolved in the water is adsorbed. The columns are removed, dried, and mounted below an electrostatic 
chamber into which the Rn from the decay of trapped Ra is continuously flowed by a stream of N 2 gas. The subsequent 
decay of Rn gives charged Po ions which are swept by the electric field onto a solid-state a counter. The content of Ra in 
the water is inferred from the measured decay rates of 212 Po, 214 Po, 216 Po, and 218 Po. The Ra extraction efficiency is >95%, 
the counting efficiency is 24% for 214 Po and 6% for 216 Po, and the method can detect a few atoms of 224 Ra per m 3 and a 
few tens of thousands of atoms of 226 Ra per m 3 . Converted to equivalent equilibrium values of the topmost elements of the 
natural radioactive chains, the detection limit in a single assay is a few times 10~ 16 g Th or U/cm 3 . The results of some 
typical assays are presented and the contributions to the systematic error are discussed. 



Keywords: radioactivity assay, water purification, solar neutrino, SNO 
PACS: 29.40.-n, 26.65.+t, 81.20.Ym 



1. Introduction 

All experiments to detect solar neutrinos have ob- 
served significantly fewer than are predicted by well 
calibrated solar models. It is now recognized, after 
publication of the first results of the SNO experi- 
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ment [1], that this deficit is the result of the transfor- 
mation of the electron neutrino into other neutrino 
species. A crucial element of the SNO experiment is 
the verification that the rates of background processes 
that might be misinterpreted as neutrino events are 
significantly less than the observed signal. We de- 
scribe here one of the procedures used to measure the 
radioactivity content of the water in the SNO detector 
and thus to prove that the background is adequately 
low. 

This article is structured as follows: after giving a 
general description of the SNO experiment, the vari- 
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ous aspects of the background measurement technique 
based on Mn-coated beads are described. Then the 
methods for measuring the Ra extraction efficiency 
and the Po counting efficiency are considered, fol- 
lowed by a discussion of systematic errors. Finally 
representative results of using the method to measure 
the Ra concentration in the heavy and light water of 
the SNO detector are given. 

1.1. Overview of the SNO detector 

The main neutrino target of the Sudbury Neutrino 
Observatory [2] is 1000 tonnes of D 2 contained in 
a transparent spherical acrylic vessel of 12 m diame- 
ter. An array of 9438 inward-looking photomultiplier 
tubes with light collectors are located on a nearly 
water-tight icosahedral support structure at ~3 m dis- 
tance from the central vessel. They see the Cherenkov 
light from neutrino interactions in the D 2 and in part 
of the 1700 m 3 of H2O which fills the volume between 
the acrylic vessel and the photomultipliers. An addi- 
tional 5700 m 3 of H 2 between the photomultipliers 
and the cavity wall provide further shielding from ex- 
ternal radiation. 

Solar neutrinos are detected in SNO through three 
distinct reactions: 

v e +d — >e~ + p+ p (CC), 
v x + d — > v x + p + n (NC), 
v x + e~ — > v x + e~ (ES), 

where x denotes any of the active neutrino species 
e,//, or t. The charged-current reaction (CC) observes 
only electron neutrinos, the neutrino type produced in 
the Sun. The neutral-current reaction (NC) detects the 
total flux of active neutrinos and has equal sensitivity 
to all flavors. The elastic-scattering reaction (ES) also 
detects all active neutrinos, but is dominantly sensitive 
to neutrinos of the electron type. 

1.2. Upper limits on radioactive contamination 

The most restrictive limits on D2O purity arise 
by considering the background for the neutral cur- 
rent reaction. The main unvetoed background for this 
reaction is the photodisintegration of the deuteron, 
d(y, p)n, whose threshold is 2.22 MeV. The only 
commonly occurring natural isotopes which emit y 
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rays above this energy are near the end of the Th 
and U decay chains, viz., 208 T1, 99.8% of whose de- 
cays produce a 2.614-MeV y ray, and 214 Bi, 2.15% 
of whose decays are accompanied by a 7 ray with 
energy greater than 2.22 MeV. If we require that the 
rate of neutron production by radioactive background 
be less than 1 per day, which is less than 10% of the 
anticipated signal, then, based on our Monte Carlo 
for the response of the detector to internal y rays, 
we set the goal of restricting the decay rate of 208 T1 
and 214 Bi in the D 2 of the SNO detector to less 
than 480/day and 32 000/day, respectively. Assuming 
radioactive equilibrium in the Th and U chains, this 
leads to upper limits of 4.2 x 10~ 15 g Th/cm 3 D 2 
and 3.3 x 10~ 14 g U/cm 3 D 2 0. 2 At these limits there 
is an average of 1.5 atoms of 208 T1 and 640 atoms of 
214 Bi in the entire D 2 volume. 

The requirements on radiochemical purity of the 
H 2 in the SNO detector are determined by the ingress 
of y rays into the D 2 region that are produced by the 
decay of 208 T1 or 214 Bi in the surrounding light wa- 
ter. Based on our early simulations for the transport of 
y rays and assuming equilibrium, this leads to upper 
limits of 3.7 x 10~ 14 g Th/cm 3 H 2 and 4.5 x 10~ 13 
g U/cm 3 H 2 for the water in the region between the 
acrylic vessel and the photomultiplier tubes, consider- 
ably less severe than required for the D 2 0. 

Four complementary techniques have been devel- 
oped by SNO to determine if the detector water meets 
these specifications. Three methods are based on flow- 
ing water to an external extraction system where a 
chemical separation of Ra, Th, or Rn is performed; the 
other method uses the Cherenkov light signals from 
the photomultipliers to directly infer the concentra- 
tions of 214 Bi and 208 T1. The first of the extraction 
methods, the MnO x assay method, is the subject of 
this paper. The second radiochemical method, which 
is described in a companion paper [3], flows water 
over a filter coated with an adsorber of hydrous tita- 
nium oxide. After water flow, the filter is eluted with 
acid to remove the extracted Ra, Th, and Pb. These el- 
ements are then concentrated, mixed with liquid scin- 



2 It is of course not true that the Th and U chains are in equilibrium 
in the water of the SNO detector all the way back to long- 
lived 232 Th and 238 U, but it is conventional in low-background 
counting to express radioactivity measurements as if equilibrium 
were present, and we follow this practice here. 



MnO K assay method 



14 January 2003 



MnO x assay method 



3 



tillator, and the /3-a coincidences of 212 Bi- 212 Po and 
214 Bi- 214 Po are detected with a photomultiplier. In the 
third chemical assay method, water from the detector 
is flowed through a degasser to liberate Rn. The Rn is 
purified and collected and its alpha decays are counted 
in a Lucas cell scintillator chamber on a photomul- 
tiplier [4]. This method can detect only the U-chain 
isotope 



222 Rn. 



2. The MnO x assay method 

In overview, water is passed through columns that 
contain beads coated with a manganese oxide com- 
pound. The coating extracts Ra from the flowing wa- 
ter, and, to a lesser extent, other dissolved species, 
such as Th, Pb, etc. [5]. After a large volume of wa- 
ter has passed through the columns, they are removed 
and dried. The Rn produced from Ra decay is swept 
from the columns into an electrostatic chamber where 
it decays. The charged Po ions from the decay of Rn 
are carried by an electric field onto an alpha counter 
where the decays of the Po are detected. 

In this section we describe all the components of the 
MnO x assay method as used in SNO, from the MnO x 
coating onto acrylic beads, the column which holds 
the beads during an assay, the system that flows water 
through the column, the apparatus for counting the 
column, the physical interpretation of the data, and the 
way in which the data are analyzed. 

2.1. Beads and Mn coating 

Solid spherical acrylic beads of 600 fim nominal 
diameter were chosen as the support material. These 
beads have good mechanical strength, little water up- 
take, low intrinsic content of radioactive elements [6], 
good resistance to the production of fine particulate 
material (fines), and good Rn emanation efficiency. 

Tests of Ra extraction efficiency and retention, Rn 
emanation, and fines production were carried out with 
two forms of coatings, designated MnC>2 and MnO x . 
The MnC>2 coating was produced by prolonged oxi- 
dation of the beads with sodium permanganate in the 
presence of sulfuric acid at high temperature. This 
yields a black coating on the bead surface. The MnO x 
coating is obtained by incomplete oxidation with the 



same reagents, and gives a dark brown coating. It was 
selected for our purpose as the Rn emanation efficiency 
was higher than with MnC>2. The exact stoichiometry 
has not been determined so we call it MnO x . 

The procedure for bead production with MnO x coat- 
ing is described in [7]. Electron microscopy reveals 
that the coating consists of a base layer that is paral- 
lel to the surface and needles of (0.3-0.5) /an length 
oriented approximately perpendicular to the surface. 
The average coating density is 0.45 + 0.1 g Mn/m 2 . 
The radioactive background of one sample of coated 
beads was measured with a Ge y-ray detector [8] to 
be <90 ng Th/g beads and 29 ± 14 ng U/g beads. 

The ability of the beads to retain adsorbed Ra was 
measured by spiking a column with 226 Ra, measuring 
its initial activity A(0), flowing a volume of water V 
through the column, and measuring the activity A(V) 
that remained on the column. These tests were con- 
ducted on small-scale analogues of our columns used 
in extraction with flow rates of 20 column-volumes 
per min and the activity was measured with a Ge y- 
ray detector. The results of these measurements can 
be described in terms of an effective distribution co- 
efficient Ki defined as 



K d (V) = 



v/v,, 



eads 



A(0)/A(V) - 1 ' 



(1) 



where Vb ea ds is the volume of the beads in the column. 
Up to V/ V b e a d s = 170 000, K d was found to be >10 6 
with an error of approximately 25%. This implies there 
is little concern regarding loss of Ra from a one liter 
column at a flow rate of 20 {/min 3 up to water volumes 
of at least 170 k(. 

Not only do the beads adsorb Ra, but they also have 
an affinity for other elements, such as Ba. By flowing 
a large volume of Ba solution over a fixed volume of 
MnO x -coated beads and measuring the Ba concentra- 
tion of the eluate, the capacity of the coated beads to 
hold Ba was determined to be ~300 mg Ba/£ beads. 
The capacity for Ra should be comparable as it also 
is a divalent atom and has only slightly greater radius. 
The other divalent atoms and monovalent atoms have 
a smaller affinity for MnO x . The water in the SNO de- 
tector does not contain such large quantities of impu- 
rities that there is any problem with the bead capacity. 



3 To avoid confusion with the number 1, we use the symbol "(" 
as an abbreviation for 'liter.' 
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2.2. Column 

During a water assay and subsequent counting, 
the coated beads are held within a column made 
from polypropylene. The internal column volume is 
-950 cm 3 and it is filled with 700.0 ± 0.2 g of coated 
beads which contain approximately 3.5 g of Mn. 

Before a column is used it is dried and counted for 
background. This is done in the same way as after the 
column has been used in an experiment and will be 
described below. We call this measurement the column 
'blank.' When counting is finished, which takes at least 
10 days, the column is filled with nitrogen, its ends are 
capped, and it is brought underground for the assay. 

To maintain the isotopic purity of the heavy wa- 
ter, columns for a D2O assay are deuterated before 
use. This is done by flowing ~3 I of D2O through the 
column with a peristaltic pump at 150 cm 3 /min un- 
til the density of the water that exits the column is 
>1.103 g/cm 3 . (The density of pure D 2 at 20 °C is 
1.105 g/cm 3 .) To reduce the entry of fines into the wa- 
ter, columns are then rinsed at a rate of 20 {/min for 
15 min. 

2.3. Wafer flow 

Most H2O assays use a sampling point halfway be- 
tween the acrylic vessel and the photomultipliers as 
it is this region of H2O that contributes most to the 
neutral-current background. Assays of the D2O usu- 
ally draw water from the vessel bottom and return it 
near the vessel top. The pH of the H2O is 5.8-6.0 and 
the pD of the D2O is 5.5-5.8. All water entering the 
detector is cooled; the temperature of the D2O and the 
immediately surrounding H 2 is (10.5-11) °C. 

After the column is attached to the water system it 
is evacuated so that the ingress of 222 Rn to the sys- 
tem is minimized. The usual water flow rate through 
an MnO x column is 20 {/min. To allow the sampling 
of very large volumes of water, and thus to achieve 
greater sensitivity, assays of the D 2 are usually made 
with 4 columns through which the water flows in par- 
allel with 20 f/min through each. 

All water that has flowed through an MnO x column 
is passed through ultrafilters that remove tiny partic- 
ulates (fines). If MnO x fines were to be carried into 
the acrylic vessel, they would bind there with Ra and 



remove it from future assays, thus leading to an under- 
estimation of the true Ra content of the water. The ul- 
trafilration units have a 3 kiloDalton molecular weight 
cutoff (approximately 3 nm). Their permeate flow re- 
turns to the water system; their concentrate flow is 
passed through a 0. 1 fim filter and then returned to the 
feed stream of the MnO x column. 

Since a very large volume of water is passed through 
the column, there is a concern that some Mn may be 
dissolved and enter the water. An indication that some 
loss of Mn occurs is provided by the visual observa- 
tion that the bead color sometimes changes from its 
initial dark brown to light brown after water flow. This 
effect is most evident on the end of the column where 
the water enters. To address this question, experiments 
have been made in which H2O at pH 7 was recirculated 
from aim 3 tank through a column and the concen- 
tration of Mn measured as a function of flow volume 
by inductively-coupled plasma mass spectrometry [9]. 
Very little Mn was found, even up to water volumes 
of 115 m 3 . Using this data and measurements of the 
mass of Mn remaining on the beads after their use in 
an extraction, it appears that a few tens to hundreds 
of mg of Mn are removed from a column in a large 
volume assay. Most of this Mn is in the form of fines 
which are blocked from entering the detector by the 
ultrafilter that follows the MnO x columns. Very little 
of the released Mn is in dissolved form: measurements 
of the D2O in the acrylic vessel show an increase in 
Mn concentration from an initial value of 0.5 ng/cm 3 
to 2.0 ng/cm 3 during the pure D2O phase, a period in 
which 61 MnO x columns were used in assays. This 
increase in Mn concentration has not deleteriously af- 
fected the operation of the SNO detector. 

2.4. Column drying and counting 

When an assay has been completed, the columns are 
removed, their ends are capped, and they are brought to 
a surface laboratory where they are dried and counted. 
Either air or boil-off gas from liquid N2 is used for 
drying. The gas is filtered and heated to 60 °C. A vol- 
ume of ~25 m 3 is used to dry a wet column from an 
assay; ~10 m 3 is sufficient for an initially dry column, 
as when counting a blank before use in an assay. 

The Ra activity on the columns is measured by at- 
taching the dried column to the gas flow loop on an 
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69) volts. The diode is held at a potential of 1000 volts 
below that of the aluminum chamber. The output of 
the diode is amplified, digitized in 1024 channels, and 
written to disk every 3 hours. Counting is continued 
for (10-30) days. Long counting times are essential in 
D2O assays to well determine the 228 Th activity. 

Fig. 1 shows the counting configuration for a single- 
column assay. When 4 columns are used they are at- 
tached to the gas flow loop in parallel in the same po- 
sition as with one column. The gas is circulated with a 
diaphragm pump at a typical flow rate of (0.3-0.4) STP 
f/min. The filter just below the ESC blocks the entry 
of fines to the chamber and is the major impediment 
to gas flow. 

Most of the Ra from a water assay is located near 
the top of the column where the water enters. To give 
the highest efficiency for the transport of Rn to the 
active ESC volume, especially the short-lived 220 Rn, 
the gas flow through the column is thus opposite to 
the previous direction of water flow. 

Energy spectra from the a detector at the end of 
counting of a typical assay are given in Fig. 2. The 
lower panel shows the spectrum of the blank before 
passing water through the column; the upper panel 
shows the spectrum of the sample after water flow. Be- 
cause the a particles lose energy as they pass through 
the surface dead layer of the diode, each Po peak 
falls gradually on the low-energy side, but much more 
steeply on the high-energy side. 



Recirculation pump 



Fig. 1. Schematic diagram of gas flow through MnO x column and 
ESC during counting. 



2.5. Interpretation of energy spectra 



electrostatic counter (ESC), as shown in Fig. 1. To re- 
move Rn and residual water vapor, the counting loop 
is initially filled with dry N2 from a compressed gas 
bottle and evacuated several times. The loop is then 
pressurized to 26 + 1 mbar with dry N2 which acts as a 
carrier gas to transport Rn from the column to the ESC. 
The total elapsed time between the end of extraction 
and the start of counting is typically (10-12) hours. 

The internal shape of the ESC chamber was de- 
signed [10] to sweep the positively charged Po ions 
from Rn decay onto an alpha detector at the apex of the 
chamber with high efficiency throughout its volume. 
The alpha counter is an 18 mm by 18 mm windowless 
silicon photodiode operated at a bias voltage of (60- 



Since it is Po produced by the decay of Rn that is 
detected, we look backward from Rn along the Th and 
U chains to the last isotope whose half life is con- 
siderably longer than the usual (10-15) day counting 
period. This leads to 1.9-y 228 Th in the Th chain and 
1600-y 226 Ra in the U chain. If any prior elements in 
either chain are present, their influence is effectively 
blocked by these long-lived isotopes. 

The relevant part of the Th decay chain and the en- 
ergies of the a particles that follow Rn is thus 228 Th 



1.91 y 



224 



Ra 



220 



3.66 d 



Rn 



216 



212 



Bi 



212 



Po 



55.6 s 
8.8 MeV a 



Po 



6.8 MeV a 



212 



0.15 



Pb 



208 



Pb, which is sta- 



10.6 h 1.01 h 298 ns 

ble. The next to last decay here takes place 64% of the 
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ble. Fig. 2 shows the major peak regions due to 218 Po 
and 214 Po in a U-dominated spectrum. Although the 
210 Po peak is also a part of the U chain, it cannot be 
used in analysis as it represents the gradual accumula- 
tion on the diode of long-lived 210 Pb from the counting 
of previous assays, tests, etc. 

Because the concentration of 226 Ra as an impurity 
in the ESC loop is higher than that of 228 Th, and the 
counting efficiency for 222 Rn is higher than that of 
220 Rn, all low-level spectra are U-dominated, rather 
than Th-dominated. Thus, the relative heights of the 
peaks in a single spectrum do not reflect the relative 
concentrations of 226 Ra and 224 Ra in the water. 

Note in Fig. 2 that the 216 Po and 2I2 Po peaks are 
somewhat higher in the sample spectrum than in the 
blank spectrum. This increase occurs because the col- 
umn extracted 224 Ra from the flowing water. In con- 
trast, the 218 Po and 2I4 Po peaks are higher in the blank 
than in the sample. This difference is not because sig- 
nificant 226 Ra was lost during water flow; rather, it is 
due to the presence of more 222 Rn in the ESC loop 
when counting of the blank was initiated than when 
counting of the sample began. 



Fig. 2. Typical energy spectra from a detector at the end of 
counting of a 4-column D2O assay. Lower panel: before assay, 
the 'blank.' Upper panel: after passing 355 k£ of D2O through 
columns, the 'sample.' The peaks in order of increasing energy 



^ 2I0 Po at 5.3 MeV, 218 Po at 6.0 MeV, 216 Po at 6.8 MeV, 214 F„ 
at 7.7 MeV, and 212 Po at 8.8 MeV. The vertical lines show the 
energy windows assigned for each peak. 



time; the other 36% proceed through 212 Bi 6,1 MeV 
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Pb. The major peaks from Th are thus 

216 Po and 212 Po, as labeled in Fig. 2. The 6.1-MeV a 
from 212 Bi makes an unresolved peak near the upper 
edge of the much stronger 21 8 Po peak. The window for 
212 Po decay is set very wide to include the detection 
of the 2.2-MeV endpoint energy electron from 212 Bi 
decay in coincidence with the 8.8-MeV a from 212 Po 
decay. 

The section of the U chain of interest is 226 Ra 
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2.6. Data analysis 

The counting data are analyzed in three steps: In 
the first step, the energy spectra are fit, energy win- 
dows are set about the four Po peaks, and the number 
of counts in each energy window is determined as a 
function of time. The second step is to fit these time 
spectra to determine the decay rates of the isotopes 
that produce the counts. These two steps are done sep- 
arately for the blank and for the sample. In the third 
step, the isotope decay rates for sample and blank and 
the water flow rate are used to determine the concen- 
tration of 224 Ra and 226 Ra in the water. 

2.6.1. Energy spectrum analysis 

The final accumulated spectrum is first passed 
through a peak recognition algorithm which identifies 
the various peaks. Once the peaks are recognized, a 
fit is made to the two principal peaks, either 212 Po 
and 216 Po for a Th-dominated spectrum or 214 Po and 
2I8 Po for a U-dominated spectrum. The energy scale 
is then set by a linear interpolation from the two 
principal peak positions and the four energy windows 
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Table 1 



Data, fitted, and calculated parameters 


for spectra in 


Fig. 2. 


Value 


Item 


Blank 




Live time (days) 


11.522 


18.288 


Counts in 218 Po window 


485 


683 


("'mint*; in Pr» winHnw 


45 


122 


Counts in 2 ^Po window 


475 


641 


Counts in 2l2 Po window 


32 


93 


Peak height H (counts) 

218p G 


49.3+j' 
63.0+ 4: ^ 


61.6+" 


214p G 

Peak location C (channels) 


74 9+ 4: ] 

' ' -4.4 


218p 


397.4+03 


396.1+" 


214p 

Width of Lorentzian W L (keV) 


513.0+o| 


512.3+jj| 


below 218 Po 
below 214 Po 
Width of Gaussian W G (keV) 


134 7+ 11 - 2 

-10.3 

105.8+J 4 


145.1+'°/ 
123.5+| 4 


above 218 Po 


84J -68 


102.2+f' 


above 214 Po 
Overlap of 214 Po into 216 Po (%) 


57.0+fl 
1.4 


60.2+1° 
1.7 


Window efficiency for 218 Po (%) 


95.5 


95.0 


Window efficiency for 216 Po (%) 


96.2 


95.4 


Window efficiency for 214 Po (%) 


97.2 


96.6 


Window efficiency for 212 Po (%) 


98.0 


98.0 



are set based on the fitted peak parameters. Finally, 
the number of counts in each window in each 3-hour 
interval is written to a file for subsequent analysis. 

This process is illustrated with the sample energy 
spectrum in Fig. 2. The principal peaks are 214 Po and 
218 Po, so this is a U-dominated spectrum. The solid 
line is the fit to the 214 Po and 218 Po peaks. Table 1 gives 
the fitted parameters and other data for these spectra. 
The 218 Po window is set so that its lower limit removes 
the 210 Po peak and its upper limit fully includes the 
212 Bi contribution. 

This figure also shows one practical complication: 
although the a resolution is good, it is not high enough 
to totally separate the peaks. As a result, each peak 
makes a contribution to the peaks at lower energy [10]. 
The extension of the fitted line in Fig. 2 below the 
214 Po peak indicates the overlap of 214 Po into the 216 Po 
energy window. Usually (1-2)% of the counts in each 
window appear as overlap counts in the adjacent lower 
energy window. Although this fraction is small, over- 
lap can be an appreciable effect. For example, since 
the 214 Po peak of the spectrum of the blank in Fig. 2 
contains 475 counts and the overlap is 1.1%, approx- 
imately 5 of the counts in the 216 Po window are from 
overlap, more than 10% of the total of 45 counts in 



this window. 

A fit function that was found to adequately approx- 
imate the shape of each Po peak is a Gaussian above 
the maximum and a Lorentzian below the maximum. 
If the peak height is H and the center (maximum) is 
at C, the number of counts m at energy E is then 

\2n 



m(E) = H 



exp 



Ue-cY 

(WJ2) 2 



for E > C 



for E < C, 



(2) 



I (E - C) 2 + (W L /2) 2 

where Wc is the width (standard deviation) of the 
Gaussian and W L is the width (FWHM) of the 
Lorentzian. Based on the Poisson probability (m\l) 
for detecting m counts in an energy bin that contains 
/ counts, (m\l) = e~ m m! //!, the likelihood function X 
is defined as 



-m(E,) [ 



l(Ei) l(Ei> 

KEiV. ' 



(3) 



where the product is over all the energy bins i in the 
window. In practice the fit is made by searching for the 
set of parameters that minimize the negative logarithm 
of £, 



■ln£ = 2[m(£,) - l(Ei)\nm(Ei)l 



(4) 



where the constant term ln[/(£,-)!] has been neglected. 
The two peaks are fit independently and three passes 
are made through the fit routine, with each pass grad- 
ually refining the four peak parameters H, C, Wg, and 
Wl, and resetting the window limits. After the last 
pass, 68% confidence regions on the parameters are 
calculated by finding the values of each parameter that 
increase - In X by 0.5, all other parameters being max- 
imized. 

Finally, the number of counts in each of the 4 energy 
windows for each time interval is determined, the best 
fit parameters are used to calculate the overlap and 
energy efficiency for each window, and all of these 
calculated values are written to an output file. 

2.6.2. Time spectrum analysis 

Examples of time spectra from the last step of en- 
ergy analysis of the data in Fig. 2 are shown in Fig. 3 
for the 226 Ra-originated components and in Fig. 4 for 
the 228 Th-originated components. Fits to these spectra 
are given by solid curves with dashed lines above and 
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Fig. 3. Time spectra of 214 Po (left side) and 218 Po (right side) for the same 4-column D2O assay as in Fig. 2. Lower panels: blank before 
assay. Upper panels: sample after passing D2O through columns. The data points are indicated by horizontal lines of 1-day duration. Since 
the count rate is very low, the data from every eight 3-hour data collection intervals have been combined. The thick solid line is the fit to 
the raw data with the encompassing dashed band indicating the approximate 68% confidence range. The major components of the fit are 
shown by the thin solid lines. See text for further explanation. 



below the fit that indicate the range of the fit based on 
the statistical 68% confidence limits for the derived 
parameters. The fit to the spectra contains four compo- 
nents for 214 Po and 216 Po, two components for 212 Po, 
and six components for 218 Po. The major components 
of the fit are separately shown by the rising and falling 
thin solid lines. The physical basis for these fits and the 
time functions used in fitting [11] will now be briefly 
described. 

2.6.2.1. U chain The U chain is the simpler case and 
will be discussed first. The observed Po isotopes, first 
218 Po and then 214 Po, are produced by two sources: 
• The decay of 222 Rn that was entrapped during expo- 
sure of the column and sections of the ESC counting 
loop to air. Neglecting any permeation of Rn into 



the loop, the decay rate of this component begins at 
some positive value and falls with the 3.82-day half 
life of 222 Rn. At time t the instantaneous decay rate 
has the form «222« ■~ Alllt , where 0222 is related to the 
initial number of 222 Rn atoms and X222 is the decay 
constant of 222 Rn. This component is shown in the 
spectrum of 218 Po by the thin falling line in Fig. 3. 
The initial amplitude 0222 depends on the extent of 
pumpout before the start of counting. Note that in 
this example there was more trapped Rn in the blank 
than in the sample. The time development of this 
component is somewhat different in the spectrum 
of 214 Po-there is a delay at the start of counting be- 
cause of the need to build up the daughter elements 
26.8-min 214 Pb and 19.9-min 214 Bi. This leads to an 
initial rise in 214 Po count rate before the fall begins. 
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Fig. 4. Time spectra of 212 Po (left side) and 2l6 Po (right side) for the same 4-column D2O assay as in Fig. 2. See caption for Fig. 3 for 
further explanation. 



• The decay of 226 Ra. 226 Ra originates from two 
sources: (1) Ra extracted from the water and (2) 
contamination with Ra of all the components of 
the ESC loop, viz., the beads, the MnO x coating, 
the column, the detector walls, etc. These sources 
are separated by counting the column both before 
and after extraction and taking the difference. This 
component has the time dependence in the 218 Po 
spectrum «226(1 - e~ Al22 ') where a 226 is related to 
the initial number of 226 Ra atoms. It is shown by 
the thin rising lines in Fig. 3. 



2.6.2.2. Th chain The parent isotopes that produce 
the observed 216 Po and 212 Po counts are 228 Th and 
224 Ra. These spectra also originate from two sources: 
• The decay of 224 Ra that (1) was produced by the 
decay of 228 Th and (2) was extracted from the water. 
Since it is the latter that we wish to determine, the 
column is counted before and after extraction. In 



the spectrum of 216 Po this component has the time 
dependence a224<?~ i224 ' where 0224 is related to the 
initial number of 224 Ra atoms and ^224 is the 224 Ra 
decay constant. It is shown in Fig. 4 in the spectrum 
of 216 Po by the thin falling line. In the spectrum of 
212 Po it is the component that gradually rises to a 
peak at about 1.5 days after the start of counting 
and then falls with the 224 Ra half life. The initial 
rise in the 212 Po spectrum is because of the delay in 
production of the 10.6-h daughter 212 Pb from 224 Ra. 
• The decay of 228 Th. The 228 Th is present as a con- 
taminant in the constituents of the ESC loop and, 
to a small extent, may be extracted from the water. 
This term is shown in the spectra of Fig. 4 by the 
thin rising line. In the 216 Po spectrum the time de- 
pendence is (3228(1 - e~ Alu ') where 0228 is related to 
the initial number of 228 Th atoms. The time depen- 
dence in the 212 Po spectrum is again delayed by the 
need to build up 212 Pb. 
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2.6.2.3. Estimation of activities 

Both decay chains are examples of the general case 
N\ — > N2 — > N3 — > ...Ni — > ... whose solution 

ii A2 /I3 A, 

is [12] 

i,M(f) = Ni (0) J] e-^Aj f] j4r ' < 5 > 

7=1 k=i k 1 

k*j 

where Ni(t) is the number of atoms of species i at time 
t and the initial condition yV,-(0) = for all i ± 1 is 
imposed. The number of decays <f, of species i during 
a counting interval that begins at time ?b and ends at 
time f E is thus 

A i N i dt = N 1 (0)}] A j pi p ( 6 ) 

where Aj = e~ A '' B - e~ x > tlL and P'. is the product in 
Eq. (5). Introducing the detection efficiency e, for iso- 
tope i and defining the activity A,- of isotope i to be 
A; = AiNi, the number of detected counts ra,- of species 
;' due to parent isotope 1 during the counting interval 
from f B to r E is 

m i (t B ,t E ) = s i ^-f j A j P i j . (7) 

Defining A222 as the initial activity of 222 Rn and 
A226 as the initial activity of 226 Ra, Eq. (7) is applied 
in the U chain to obtain the number of counts of 218 Po 
and 214 Po in each counting interval. Similarly, defining 
A 224 as the initial activity of 224 Ra and A228 as the 
initial activity of 228 Th, Eq. (7) is used in the Th chain 
to obtain the number of counts of 216 Po and 212 Po 
in each counting interval. Because their half lives are 
much shorter than the 3 -hour interval on which data is 
recorded, we make the approximation that the 55.6-sec 
22() Rn and 0.15-sec 216 Po decays occur immediately 
after the 224 Ra decay. 

In addition to these major terms, there are several 
additional components, whose amplitudes are rela- 
tively small, which also contribute to the time spectra: 
Except for 212 Po, each Po isotope has an overlap con- 
tribution from the adjacent Po isotope at higher en- 
ergy in the spectrum. These terms are added onto the 
expression for m,(f B ,fE), using for 216 Po the overlap 
fraction measured when Rn gas is added to the ESC 
loop, and for 214 Po and 218 Po the overlap fraction mea- 
sured when a Th source was in the ESC loop. Finally, 



for 218 Po, a term is added to m, to take account of the 
contamination of 212 Bi. The relative amplitude of this 
component is approximately determined by the 36% 
branching fraction of 212 Bi to 2()8 T1 and was checked 
by fits to data with a Th source. All these additional 
terms are included in the fitted total spectra shown in 
Figures 3 and 4, but for clarity only the major terms 
are illustrated. 

For each of the four Po isotopes a likelihood func- 
tion is formed which is identical to Eq. (3) with the 
factor m(Ei) replaced by the expression for the total 
number of predicted counts m, in Eq. (7) and the prod- 
uct taken over all time intervals i. The grand likelihood 
function is defined as the product of these individual 
likelihood functions for 212 Po, 214 Po, 216 Po, and 218 Po. 
The fit is made by solving for the four variables A 222, 
A224, A226, and A228 that maximize the grand likeli- 
hood function and the parameter confidence regions 
are found using the same prescription as described 
earlier for the energy spectrum. In the spectra from a 
blank the 224 Ra and 228 Th are in equilibrium so the ad- 
ditional constraint A224 = A228 is imposed. In making 
the fits, negative values of the parameters are not al- 
lowed, since they correspond to an unphysical regime. 

The final step of time spectrum analysis is to de- 
termine the goodness of fit between the fitted spectra 
and the data. For the set of 4 spectra, the Poisson like- 
lihood x 1 goodness of fit parameter defined in [13], 
X 2 = 2 Yji\. m i ~ h + k In^iMi)] is calculated, where the 
sum is over all time intervals i and Z, is the number of 
observed counts in interval ;'. Since there are usually 
very few counts in each of the 3-hour data intervals, 
the probability distribution of this parameter is not the 
same as that of the standard x 1 distribution, but must 
be obtained by simulation. This is done by making a 
large number of simulations, each of them based on 
the best fit parameters. Each simulation is then fit to 
obtain new parameters and the value of x 1 is calcu- 
lated. This determines the distribution of x 2 for data 
with the same average number of counts per interval 
as in the real data. The probability level is then set by 
the placement of the value of x 2 from the real data 
in the distribution from simulation, a procedure anal- 
ogous to that described in [14]. 

The results of fitting the spectra in Fig. 3 and 4 are 
given in Table 2. The values used for counting effi- 
ciency are discussed in Section 4. There is consider- 
ably more 224 Ra in the sample than in the blank, im- 
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Table 2 

Data values and fitted parameters for time spectra in Figs. 3 and 4. 







Value 


Item 


Blank 


Sample 


Number of counting intervals 


92 


146 


Counting efficiency (%) 
212p 


4.4 


4.4 


214 Po 


16.4 


16.3 


216 Po 


4.1 


4.1 


218p 


15.1 


15.1 


Activity at start of counting (decays/day) 

224 Ra ' 76+'° 
228 Th 76 +10 

226 Ra 313 + \1 
222 Rn 184 +f0 
Goodness of fit probability (%) 30 ± 1.5 


278 +^ 
93 +15 

280 
69 

49 ± 1.6 



plying that 224 Ra was extracted from the water. Within 
errors the 228 Th activity is the same in both sample 
and blank, indicating that little Th was extracted. The 
disparity in 222 Rn content shows that there was much 
less Rn in the column for the sample than in the blank. 
Surprisingly, there is also slightly less 226 Ra in the 
sample than in the blank, implying that some 226 Ra 
may have been leached from the column. A loss of Ra 
of this magnitude is consistent with Eq. (1). The fits 
to both sample and blank spectra are quite acceptable, 
as indicated by the goodness of fit probabilities. 

The results of this example are very much like those 
of an average assay. Typical values for the activities 
when 4 blank columns are counted are 75 decays of 
224 Ra per day and 250 decays of 226 Ra per day, with a 
variation of 20% and 40%, respectively. These should 
be compared with measured activities when the ESC 
loop is counted with no columns present, but the loop 
closed with a dummy insert, in which case the typical 
activities are 10 decays of 224 Ra per day and (150- 
200) decays of 226 Ra per day. 

2.6.3. Concentration calculation 

We will now relate the deduced activities of 224 Ra 
and 226 Ra, A224 and A226, respectively, to the concen- 
tration CR d of Ra in the water. The first step is to de- 
termine the net activity of the sample after subtraction 
of the blank activity, defined as 

^4 net _ ^sample _ j^Uaak^^^fi^"-^) ^ 



(8) 



where ^^ d ™ ple j s the activity at the start of counting 
(SOC) of the sample and A^™ k is the activity at the 
start of counting of the blank, as calculated by the fit 



procedure in Section 2.6.2.3. The exponential factor in 
this equation decays the blank activity from the time of 
counting the blank to the time of counting the sample, 
which assumes that the blank activity is not supported 
by elements higher in the chain. This equation applies 
for both 224 Ra, in which case /Isupponing is the decay 
constant of 228 Th, and 226 Ra, in which case /Isupponing is 
the decay constant of 226 Ra. Since the time from blank 
counting to sample counting is usually much less than 
the half life of the supporting isotope, the exponential 
factor is close to unity. 

A special case arises for 224 Ra: if the 228 Th activity 
of the sample exceeds that of the blank, then there is 
evidence that 228 Th was extracted from the water. In 
this case it is necessary to subtract from A™' given by 
Eq. (8) an additional term. See [7] for details. 

During the extraction process, the number of atoms 
of Ra N extracted from the water and adsorbed on the 
column is governed by the differential equation 

dN 

— = £extractionC Ra F(f) - AN(t), (9) 

df 

where ^extraction is the extraction efficiency and F (/) is 
the flow rate. The concentration is assumed here to be 
a constant. Since our assays are made in recirculation, 
to satisfy this assumption we must return the water to 
a point in the vessel that is far from the sampling point 
and the volume of water that is sampled must be less 
than the total vessel volume. The general solution of 
Eq. (9) under the initial condition N(0) = is 



N(t) 



Jo 



F(x) dx. 



(10) 



The extraction process is often interrupted and the 
flow rate may not be constant. To model these changes, 
we break up the extraction time into / ex tr intervals dur- 
ing each of which the flow rate is constant at the value 
ft. If we define the beginning and ending times of flow 
interval i to be fbegin, and f en d,, respectively, then the 
number of Ra atoms that were extracted from the wa- 
ter and present on the column at the start of counting 

fS0C is Nfeoc) = SextractionC Ra <I>, where 



i=\ 



'^(*SOC - *eDdp _ g -/ *( f SOC -f begin ( ) 



(11) 



is called the flow-saturation factor. O has units 
of volume. From the definition of activity, A net = 
/WVfeoc). we obtain the Ra concentration 



12 



MnO x assay method 



14 January 2003 



Table 3 

Factors needed to compute concentration for example data in 
Figs. 3 and 4. 

Value 



Item 


224 Ra 


226 Ra 


Exponential factor in Eq. (8) 


0.981 


1.000 


Net activity A net (decays/day) 


197 


< 20.0 


Flow-saturation factor <D (kt) 


213.8 


320.4 


Decay constant /tR a (l(T 3 days~') 


189 


0.00119 


Extraction efficiency extraction (%) 


95 


95 


Ra concentration Cr 2 (atoms/k£) 


5.1 


< 55200 


Parent cone. (10 g Th or U/cm ) 


2.8 


< 0.06 



A net 



(12) 



£extraction^Ra c t ) 

As derived here the units of Cr 3 are Ra atoms per 
unit volume of water. It is more customary to express 
the results in terms of the parent isotopes at the top of 
the decay chains, assuming equilibrium. This implies 
in the Th chain /I224C224 ~ ^232(^232, where C232 is the 
concentration of 232 Th, and in the U chain /I226C226 = 
^238^238 where C238 is the concentration of 238 U. Fur- 
ther, we convert to the more common units of mass 
per unit volume and remove the isotopic fraction of 
232 Th and 238 U. This gives the final concentration of 
Th, Cjh, and of U, Crj, 

gTh 

(13) 



C Th = 2.851 x 1(T 15 



a net 
^224 



£extraction < l ) 224 cm J 
a net 
/i 



1 226 



gu 



Cu = 0.937 x 1(T 15 m -. 

e extraction*A'226 cm 

where the activity A is in decays per day and <5 is 
measured in k£. 

For the example set of data given in Table 2, the 
factors needed to compute Cr 3 and the results are listed 
in Table 3. Since the 226 Ra activity of the sample was 
less than the blank, we can only derive an upper limit 
for the 226 Ra activity, which we set by the quadratic 
combination of the errors for the fit of the sample and 
blank. 



3. Measurement of extraction efficiency 

In this section we describe the ways in which the 
efficiency of extraction of Ra from water by MnO x - 
coated beads has been measured. The efficiency is usu- 
ally very high, 95% or more. 

The extraction efficiency for Ra is defined as the 
probability that an atom of Ra dissolved in the water 



which flows over the column will be adsorbed on the 
column and held there until extraction ends. As dis- 
cussed below in Section 6, not all Ra that is present 
in the system may be dissolved; further there is a con- 
cern that Ra captured on the column at the start of 
water flow may be leached from the column by the 
time flow has ended, especially in assays that sample 
a very large volume of water. 

The efficiency of most extractions in the H2O is 
determined by passing the water over two columns 
(or with multicolumn assays, two sets of columns) in 
series, called here the 'upstream' and 'downstream' 
columns. Assuming both columns have the same ef- 
ficiency, which is valid because we are considerably 
below capacity and the distribution coefficient is inde- 
pendent of concentration [15], and no leaching occurs, 
then it is easily shown that this efficiency is given by 



^extraction — 1 



A net 

downstream 
A net 
upstream 



(14) 



where A net is the measured column activity after blank 
subtraction and correction for Th extraction. To help 
to satisfy the first assumption, large quantities of the 
beads, sufficient for filling many columns, are mixed, 
and the columns are always filled with the same bead 
mass. The validity of the second assumption, that no 
leaching occurs, has been tested by spiking beads 
with 226 Ra activity and flowing large volumes of wa- 
ter through them, as described above in Section 2.1. 
Insignificant loss due to leaching was observed up to 
equivalent volumes of 170 k{. Assuming the correct- 
ness of these assumptions, we obtain from Eq. (14) 
in most high-activity assays an efficiency in the range 
of (95-100)%. In the analysis we use the extraction 
efficiency from Eq. (14) whenever possible; in other 
cases an efficiency of 95% is assumed with a 5% sys- 
tematic uncertainty. 

The extraction efficiency has also been measured 
by adding weak 224 Ra spikes to ~1 m 3 of water. By 
comparison of the known initial 224 Ra activity with 
the activity inferred from counting the column on an 
ESC, these measurements showed that Ra was rou- 
tinely extracted from water by MnO x with >90% ef- 
ficiency [16]. 

The extraction efficiency of MnO x for Th at the flow 
rate of a standard assay is quite low (<10%) [17]. 
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Table 4 

Factors entering the counting efficiency when ESC is at standard operating conditions (26 mbar and 1000 V). 









Value and uncertainty (%) 


Factor 


Definition 


Th chain 


U chain 



^emanation Probability Rn atom escapes from MnO x coating 

and enters gas phase 
^volume Probability Rn atom is in ESC active volume 

when it decays 

^detection Probability Po atom from Rn decay in ESC active 
volume produces detected pulse in a counter 

^window Probability detected a from Po decay is in anal- 
ysis window for that Po isotope 



30 ± 7 ( 220 Rn) 75 ± 10 ( 222 Rn) 

83.3 ± 1.5 ( 220 Rn, 1 column) 81.9 ± 1.5 ( 222 Rn, 1 column) 
55.9 ± 1.5 ( 220 Rn, 4 columns) 55.4 ± 1.5 ( 222 Rn, 4 columns) 

26.4 ± 3.4 ( 212 Po) 40.7 ± 1.8 ( 214 Po) 
25.6 ± 3.3 ( 216 Po) 38.2 ± 1.7 ( 218 Po) 

97 ± 1 (typical, calculated in each assay for each Po isotope) 



4. Measurement of counting efficiency 

This section gives the results of counting efficiency 
measurement. For short-lived 220 Rn, the major loss of 
efficiency comes from decay before it escapes from 
the column; for long-lived 222 Rn, the major loss is dur- 
ing the process of detection. In one-column assays the 
approximate counting efficiency is 6% for the Po iso- 
topes that follow 220 Rn and 25% for those that follow 
222 Rn. 

The counting efficiency for each Po isotope is de- 
fined as the probability that if an atom of 226 Ra (for 
218 Po and 214 Po) or 224 Ra (for 216 Po and 212 Po) decays 
within the column on the ESC, then the resultant Po 
isotope will yield a detected count within the energy 
window chosen for that Po isotope. This efficiency 
was determined in two ways: (1) by 224 Ra spike ex- 
periments similar to those described earlier for mea- 
surement of extraction efficiency and (2) by measure- 
ments of all the separate factors that together make up 
the counting efficiency. 

The spike technique has the advantage that it mea- 
sures the entire throughput of the assay procedure 
from added Ra to detected Po. In these experiments, 
224 Ra activity was added to 1 m 3 of ultrapure H2O 
and the water was flowed through an MnO x column at 
20 £/min. The spike activity was determined by elec- 
troplating an aliquot of the same activity as added to 
the tank onto a planchet and counting it with an a 
detector of known efficiency. To separately determine 
the extraction efficiency, these experiments were made 
with upstream and downstream columns that were 
individually counted. Removing the extraction effi- 
ciency, which was consistently (98-99)%, the count- 
ing efficiency for 216 Po was measured [18] to be 6.0 + 
1 .7%, where the uncertainty is the standard deviation 



of the separate measurements. The large uncertainty is 
believed to be mainly due to differences in the MnO x 
coating on the beads, which lead to differences in em- 
anation of Rn. 

During early work to develop the MnO x method, 
much attention was devoted to understanding the var- 
ious factors that enter the counting efficiency. Follow- 
ing the trail from adsorbed Ra to detected Po, this ef- 
ficiency can be broken down into the terms defined in 
Table 4. A brief discussion of these factors and how 
they are measured follows. 

The emanation efficiency has been determined by 
placing beads spiked with 224 Ra or 228 Th in front of a 
shielded y counter and measuring the activity of 212 Pb 
as Ar or N2 gas was flowed at low pressure through the 
beads. The major purpose of these experiments was 
to test Rn emanation with different types of beads and 
methods of Mn coating. The values given in Table 4 
are for our present bead type and operating conditions. 
The difference in emanation efficiency between 220 Rn 
and 222 Rn is mainly due to the difference in their half 
lives. 

The volume efficiency can be calculated based on 
the known flow rate, volumes, and pressures in the 
ESC loop. At the typical flow rate the transfer time 
from the beads to the ESC is (1-2) sec with 1 col- 
umn and (4-10) sec with 4 columns. Assuming good 
mixing, the Rn remains in the ESC chamber for ap- 
proximately 25 sec, after which it returns to the col- 
umn where it either may be readsorbed on the Mn 
coating or may make a second pass through the loop. 
Summing the probability of decay in the ESC active 
volume for all possible passes leads to the values in 
Table 4. Again, the major cause of efficiency differ- 
ence between 220 Rn and 222 Rn is due to their half life 
difference. The uncertainty arises from lack of knowl- 
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edge of the extent of mixing and the probability of 
readsorption, which, from measurements of 222 Rn ad- 
sorption in air, is taken to be 4% with 1 column. As 
a check on these calculations, a 224 Ra-spiked column 
was counted by itself and then with 3 other columns. 
The ratio of the count rates in the two configurations 
was 0.655, which agrees well with the calculated ratio 
of 0.671 +0.022. 

By adding sources that emit known quantities of 
222 Rn or 22() Rn to an ESC loop, the efficiency prod- 
uct evoiumeSdetectionSwindow has been measured [10] to 
be 35 + 1.4% for 214 Po and 22 + 2.8% for 216 Po at 
our standard operating conditions. For the experimen- 
tal arrangement used in [10], the volume efficiency is 
calculated to be 90.2 + 1.5% for both 220 Rn and 222 Rn 
and the energy window efficiency is calculated to be 
95.4%. We divide by these two factors and give the 
resultant detection efficiency in Table 4. 

Because of its long half life, essentially all of the 
charged 218 Po ions from 222 Rn decay are collected pro- 
vided the ESC pressure is <100 mbar, as is almost 
always satisfied during counting. But, because of the 
short 150-ms half life of 216 Po, the daughter of 220 Rn, 
the efficiency for 216 Po detection depends more sensi- 
tively on the N2 pressure in the counting loop. Mea- 
sured values are given in [7]. The pressure in the ESC 
counting loop generally rises less than 1 mbar/day 
which leads to a reduction in the 216 Po detection effi- 
ciency of no more than a few percent during the course 
of a normal 10-day counting period. One might ex- 
pect that the detection efficiency for 212 Po would be 
approximately 64% of that for 216 Po, because of the 
64% branching ratio of 212 Bi to 212 Po. Measurements, 
however, show that the 212 Po efficiency is somewhat 
greater than that of 216 Po, presumably due to collec- 
tion of some of the 212 Pb ions if the 216 Po decays in 
flight before it reaches the diode (and perhaps to col- 
lection of some charged 212 Bi atoms which recoil into 
the gas when 212 Pb on a wall decays). Measurements 
with a Th foil in the ESC loop gave a 212 Po to 216 Po 
efficiency ratio of 1 .034. For similar reasons, the 214 Po 
efficiency is somewhat higher than the 218 Po efficiency. 
Measurements with 222 Rn in six different ESCs gave 
a 214 Po to 218 Po efficiency ratio of 1.065 + 0.006. 

The final efficiency factor, the window efficiency, 
is calculated by numerical integration over the chosen 
energy window using the fitted line shape parameters. 
Because the energy windows are set as wide as pos- 



sible, its value is quite high, typically 98% for 212 Po, 
97% for 214 Po and (95-96)% for 216 Po and 218 Po. 

By multiplying the first three factors listed in Ta- 
ble 4, this second method of efficiency determina- 
tion gives a counting efficiency in 1 -column assays of 
6.4 + 1.7% for 216 Po and 25.0 + 3.5% for 214 Po. The 
window efficiency is excluded here because it is cal- 
culated separately for each experiment and each en- 
ergy window. These are the standard counting efficien- 
cies that we use for data analysis. Including a typical 
95% energy window efficiency, the 216 Po efficiency is 
6.1 + 1.6% which agrees well with the efficiency mea- 
sured by adding 224 Ra spikes to water of 6.0 + 1.7%. 
In 4-column assays the counting efficiency is reduced 
by factors of 0.67 for 216 Po and 0.68 for 214 Po. 



5. Systematic uncertainties 

The various systematic effects are listed in Table 5. 
The largest entry is the 27% uncertainty in the counting 
efficiency for 224 Ra, discussed in Sec. 4. The other 
major entries are briefly considered here; the terms 
whose uncertainties are 5% or less are described in [7]. 

Table 5 

Sources of systematic uncertainty in Ra concentration. All entries 
are symmetric unless otherwise indicated. 



Syst. uncertainty (%) 
2^ 



Contributor 


D 2 


H 2 


D 2 


H 2 


Water flow 










Extraction efficiency 


5 


5 


5 


5 


Water sample volume 


5 


2 


5 


2 


Resampling correction 


<+4 





<+4 





Feedback correction 


3 


3 


3 


3 


MnO x fines recirculation 


+5 


+5 


+5 


+5 


Ra leaching 


+ 10 


+3 


+ 10 


+3 


Ra adsorption 


<+2 


<+2 


<+2 


<+2 


Ra background (see text) 










Data acquisition 










Counting efficiency 


27 


27 


14 


14 


Relative ESC counting eff. 


5 


5 


5 


5 


Ra distribution in column 


+5 


+5 


+ 1 


+ 1 


Drifts during counting 


2 


2 


1 


1 


Data analysis 










Overlap 


< 1 


< 1 


< 1 


< 1 


Relative line efficiencies 


3 


3 








Time functions 


5 


5 


5 


5 


Total, 


<+32 


<+30 


<+21 


<+18 


combined quadratically 


-29 


-29 


-17 


-17 
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In some high-volume assays a lower 226 Ra activ- 
ity is seen in the sample than in the blank, implying 
that some leaching may be occurring. If this is true, 
then 224 Ra must of course also be lost, but no decrease 
in 224 Ra activity is seen. This is because the ratio of 
concentrations between water and column is higher 
for 224 Ra than for 226 Ra, so the loss of 224 Ra due to 
leaching is less than the pickup from the water. Mea- 
surements very rarely show more than a 10% loss of 
226 Ra activity, so we assign a +10% uncertainty for 
Ra leaching in D2O extractions. This effect is smaller 
in H2O extractions as much less water is circulated. 

The 224 Ra background from 228 Th plated on the 
walls of the extraction apparatus of the D2O system 
was measured as follows: A closed loop was first made 
whose water flow bypassed the acrylic vessel. The 
D2O in this loop was flowed through an MnO x column 
to extract all Ra in the water, and the loop was left 
without flow for 4 days, the typical assay duration. 
Then any fresh Ra that had grown in was extracted by 
circulating the water in the loop over an MnO x column 
for 3 hours. After blank subtraction the net activity 
was measured to be 16+ 2 g decays of 224 Ra/day. This 
should be compared with the activity extracted in a 
typical assay, which is 200 + 45 decays of 224 Ra/day. 
The background from this source is thus of 
the typical assay activity. Adopting the upper limit 
as the uncertainty and converting to the equivalent 
concentration of 232 Th, we find that this systematic 
background is -0.6 x 10~ 15 g Th/cm 3 D 2 0. 



6. Results of water assays 

We give here the results of MnO x assays of the 
224 Ra and 226 Ra concentrations in the heavy and light 
water of the SNO experiment. 

The results of several 4-column assays of the D2O 
are given in Table 6. To infer the average concentra- 
tion of Th and U in the heavy water from these results 
requires an understanding of flow patterns, plating, 
solubility, diffusion, etc. and is the subject of a future 
article. Nonetheless, if one compares these results to 
the upper limit goals for radioactive purity given in 
the Introduction, it is apparent that the U value is con- 
siderably below the goal and the Th measurement is 
nearly a factor of 2 less than the goal. 



Table 6 

Equivalent Th and U content for a few typical 4-column MnO x 
assays of D2O inside the acrylic vessel. Errors are statistical with 
68% confidence. The systematic error is given in the 2nd and 
4th columns of Table 5. To convert from the units used here, 
g impurity/cm 3 D2O, to units of g impurity/g D2O, divide by the 
D2O density of 1.106 g/cm 3 at 13 °C, the temperature of the D2O 
when it flows through the columns. 



Assay 


Volume 




Date 


D 2 




(yy/mm/dd) 


(k£) g Th/cm 3 D 2 

r,A1 c A n..m-lS ,1+0.27n 


g U/cm 3 D 2 

A 1. ,10-15 ,1+O.0!h 



00/08/22 
00/10/30 
01/01/30 
01/04/16 
01/05/14 



297.0 
319.4 
381.3 
383.6 
434.1 



2.3xl0- 15 (1^-p 
2.5X10- 15 (1+° : |) 
2.8xl0-' 5 
2.1xl0-' 5 (l!n|) 
2.2xlQ-' 5 (lffi^) 



3.9x10" 



;8$ 

^ -0. 



■0.5 V 

0.0(0.0-9.6xl0~ n ) 
0.0(0.0-6.2x10-") 
1.9X10- 16 (l+^O) 
0.0(0.0-4.0x10-") 




•3 40 



5.0 10.0 
Time after start of counting (days) 

Fig. 5. Combined time spectra and fit for 212 Po (top) and 216 Po 
(bottom) for six 4-column D2O assays. The data from every eight 
3-hour data collection intervals have been combined. 



A useful way to visualize the signal of extracted 
224 Ra is to combine the data from several separate 
measurements. This can be done by simply adding 
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the counts in the separate time intervals after start of 
counting for each experiment, with appropriate cor- 
rections for times when the counter was not active. 
Figure 5 shows the 212 Po and 216 Po spectra obtained 
by combining in this manner six 4-column assays of 
D2O. The decay of 224 Ra is much more apparent than 
for the similar plot for a single measurement in Fig. 4. 

Table 7 

Equivalent Th and U content for a few typical 1 -column MnO x 
assays of H2O. Errors are statistical with 68% confidence. The 
systematic error is given in the 3rd and 5th columns of Table 5. 

Assay Volume 
Date H 2 

(yy/mm/dd) (kg) g Th/cm 3 H 2 Q g U/cm 3 H 2 Q 

Water between acrylic vessel and photomultipliers 



Upstream column 



00/01/18 


51.5 


1.0x10" 


■13 


(,+0.14) 


4.9x10-" (l«»p 


00/05/10 


47.0 


1.3x10" 


■13 


( 1+ 8:Af) 


2.8xl0-' 5 (l + n|) 
l.lxlO-' 5 (if") 


00/10/23 


110.2 


6.7x10- 


■14 




01/02/05 


45.6 


5.9x10" 


-14 


(1+0.08) 


0.0(0.0-2.4xl0- 16 ) 


01/05/23 


51.9 


6.1x10" 


-14 


/, +8:811 1 

*--n or- 1 


1.3X10- 15 (ltJ5) 



Downstream column 



00/01/18 5l5 3.6xl0- 15 (T5S) 0.0(0.0^1.2xl0- lb ) 

00/05/10 47.0 2.0xl0" 15 (l4f ) 0.0(0.0-5.6xlO- 16 ) 

00/10/23 110.2 2.3xl0- 15 (l+ofj) 0.0(0.0-1.4xl0- 16 ) 

01/02/05 45.6 7.7xKr 16 (l+fg) 0.0(0.0-2.7xl0" 16 ) 

01/05/23 51.9 0.0(0.0-1.2xKr 15 ) 0.0(0.0-1.9xlQ- 16 ) 

Water outside photomultipliers 
Upstream column 

99/05/10 55! 1.7X10-' 3 (lfpno) 1.5xKr 14 (1^) 

00/05/15 34.9 1.2xl0- 13 (l4™) 0.0(0.0-6.2xl0" 16 ) 

01/04/05 14.2 8.4xKr 14 (1™ : }°) 3.1xlQ-' 5 
Downstream column 

99/05/10 553 4.6xl()-' 4 (l!!!fj) 4.2xl0- 15 (FJnjj) 

00/05/15 34.9 5.8xl0" 15 0.0(0.0-6.1xl0" 16 ) 

01/04/05 14.2 1.3X10" 14 (1^) 1.9xlQ-' 4 0!^) 



Two types of Ra assays are often made in the H2O 
system. The first has a sampling point on the mid- 
plane of the acrylic vessel halfway between the ves- 
sel and the photomultipliers with the purified water 
returned to the midplane at two points, with the clos- 
est at 7.0 m distance from where the sample is taken. 
Since the Ra concentration in this region is higher than 
in the D 2 0, an adequate 224 Ra signal is obtained with 
a single MnO x column and a sample volume of only 
~50 k£. The results of several measurements of this 
type are given in Table 7. Both upstream and down- 
stream columns are usually used; applying Eq. (14) to 
224 Ra, it is seen that the extraction efficiency is in the 
range of (96-98)%. The other type of frequent H 2 



assay is of the water in the region outside the photo- 
multipliers. The sampling point is 4 m below the sur- 
face and the water is returned to either the midplane 
of the acrylic vessel or near the water surface. Some 
representative results are given in Table 7. Compar- 
ing the results for the region between the acrylic ves- 
sel and the photomultipliers to the upper limit goals 
for radioactive purity stated in the Introduction, one 
sees that the U content of the bulk H2O amply meets 
the goal and the Th content is slightly, but not signif- 
icantly, higher than the goal. 

For all the D2O assays in Table 6, and all but one 
of the H2O assays in Table 7, the 228 Th activity of 
sample and blank were equal within the statistical 68% 
confidence range. This is expected as the extraction 
efficiency for Th is low, and most Th, if present, is 
expected to be plated onto surfaces, or possibly be 
held in the water in the form of a soluble chemical 
complex. 



7. Conclusion 

The MnO x assay method has shown that the Ra ac- 
tivity in the bulk D 2 and the bulk H 2 of the SNO 
detector is sufficiently low that the rate of photodisin- 
tegration background events produced by Ra impuri- 
ties should be no more than 10% of the rate of neutral- 
current solar neutrino events predicted by the standard 
solar model. The sensitivity of the MnO x assay method 
is thus adequate for the needs of the SNO detector. 

Since the MnO x method is reliable and robust, SNO 
has used it to make many water quality measure- 
ments in addition to those in Sec. 6. These experiments 
mostly tested individual sections of the water system, 
such as the piping from the acrylic vessel to the MnO x 
columns, the concentrate of the reverse osmosis appa- 
ratus, the contents of various storage tanks, the purity 
of the D2O during detector filling, etc. A desirable 
feature of the method is that after extraction the col- 
umn is simply dried and counted; there is no need for 
any chemical processing, which may introduce back- 
ground from reagents and the processing system. 

Defining the detection limit as ~3 times the standard 
deviation of the background [19], then, as can be seen 
from the results in Table 2 for a typical assay, the ulti- 
mate sensitivity of the MnO x method in a single mea- 
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surement, as presently used, is ~5 x 10 16 g Th/cm 3 
D 2 and ~2 x 1(T 16 g U/cm 3 D 2 0. 
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